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D 

^ , Abstract: We apply gauge/gravity duality to compute l/N^ corrections to the heavy quark 

>- ] potentials of a quark-anti-quark pair (QQ) and of a quark-quark pair (QQ) immersed into the 
O 
O , 

, of supergravity modes between two string worldsheets stretching from the UV boundary of AdS 
space to the black hole horizon in the bulk and smeared over S . We find that the contributions to 
^ the QQ potential coming from the exchanges of all of the relevant modes (such as dilaton, massive 
scalar, 2-form field, and graviton) are all attractive, leading to an attractive net QQ potential. 
We show that at large separations r and/or high-temperature T the potential is of Yukawa-type, 
dominated by the graviton exchange, in agreement with earlier findings. On the other hand, at 
small-rT the QQ potential scales as ~ (1/r) ln(l/r T). In the case of QQ potential the 2-form 
contribution changes sign and becomes repulsive: however, the net QQ potential remains attractive. 
At large-rT it is dominated by the graviton exchange, while at small-rT the QQ potential becomes 
Coulomb-like. 
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1. Introduction 



Heavy quark potential in the vacuum of the A/" = 4 supersymmetric Yang-Mills (SYM) theory 
[1,2] was one of the first results obtained using the Anti-de Sitter space/Conformal Field Theory 
(AdS/CFT) correspondence [3-7]. The first attempt to generalize the calculation [1,2] to the case of 
the A/" = 4 SYM medium at finite temperature (T) was made in [8,9] shortly thereafter. The authors 
of [8, 9] identified two string configurations contributing to the finite-T heavy quark potential in 
the AdS/CFT framework, shown below in Fig. ^ and calculated their contributions. Defining the 
heavy quark potential as the free energy of the quark-anti-quark system immersed in the thermal 
bath, and assuming that the two string configuration from Fig. ^ contribute to the free energy on 
equal footing, the authors of [8,9] obtained a potential which had a kink when plotted as a function 
of the quark-anti-quark separation r. The kink (a derivative discontinuity) was due to a transition 
from the regime when one string configuration dominated the potential to the regime where another 
one was more important. 

The kink feature of the result obtained in [8, 9] appears to be an artifact of the large- A''^ 
approximation, as was noted in [10,11]. In [10,11] it was argued that there are of the order of 

configurations with two straight strings (shown in the right panel of Fig. such that the two 
configurations in Fig. |^ do not come in on equal footing, contrary to the assumption made in [8,9]. 
Here is the number of D3 branes, and configurations may be interpreted as resulting from 
Chan-Paton indices of the two strings, each string having Nc of indices (see also [12], where A"^ 
configurations emerge from the topological x Zn symmetry of thermal SYM theory). This large 
number of string configurations enhances the terms which would be otherwise subleading. As was 
pointed out in [11], while the exchanges of supergravity fields between the straight strings shown 
in Fig. 1^ are normally A|-suppressed, they become leading-order due to the fact that there are A"^ 
string configurations to enhance their contributions. The primary goal of [11] was determination of 
Debye mass, defined as the screening mass for chromo-electric modes in the plasma following [13]: 
therefore, for the heavy quark potential, the authors of [11] used the glueball masses calculated for 
QCD3 from AdS/CFT correspondence in [14,15], to obtain only the exponential form of the large-r 
asymptotics of the contribution coming from the supergravity field exchanges to the heavy quark 
potential. 

The question of the definition of the heavy quark potential at finite-T in the context of AdS/CFT 
was raised in [16]. There it was pointed out that one can define two heavy quark potentials: 
color-singlet and color-adjoint [17,18]. It was further conjectured in [16] that the hanging string 
configuration (left panel in Fig. ^ may correspond to the color-singlet potential, while the two- 
strings configuration (right panel in Fig. |^) may give the adjoint potential (which is zero in the 
large- A'c limit). It was further suggested in [16] that the singlet potential at large-r T, for which 
the solution corresponding to the hanging string configuration becomes complex, may be obtained 
by analytically continuing the solution into complex domain. Such analytic continuation led to an 
absorptive potential which had both real and imaginary parts, with the real part falling off as a 
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power of r at very large-r, in some similarity to the weak-coupling perturbative results [19-25]. The 
issues associated with different definitions of the heavy quark potential are presented below in Sec. 
0, where we present both definition of the finite-T heavy quark potential, and discuss the pros and 
cons of using each of them. 

The central goal of this paper is to find the full contribution to the heavy quark potential 
coming from the exchanges of supergravity fields between two straight strings oriented in opposite 
directions, as shown in Fig. ^ below. In the language of [16] this contribution would correspond to 
the color-adjoint potential. The calculation is presented in Sec. ^ and is carried out in Euchdean 
time, making it free of subtleties of the real-time formalism of thermal field theories. We concentrate 
on the case of the strings smeared over the 5*^, which correspond to Polyakov lines defined below in 
Eq. ( |3.1| ) averaged over the directions of their coupling to six scalars of A/" = 4 SYM. In this case 
the supergravity exchanges are limited to k = Kaluza-Klein (KK) modes. The contributions of all 
the relevant supergravity modes (dilaton, massive scalar, 2-form field and graviton) are calculated 
in detail in Sec. ^. We find that all the exchanges give attractive contributions to the QQ potential. 
The net result for the QQ potential in momentum space is shown in Fig. (lower line). The 
potential turns out to be attractive at all values of rT (which is the only relevant dimensionless 
parameter in the problem). The large-r T asymptotics of the net QQ potential is given by the 
exponential decay of Eq. (p.97|) , with the exponent determined by the first pole of the momentum 
(g)-space potential along the imaginary-g axis (in agreement with [11]). The pole corresponds to the 
glueball mass in QCD3 calculated in [14,15]. We also find the residue of the pole (the factor in front 
of the exponential). The small-rT asymptotics of the potential is rather peculiar, since it is not 
quite Coulomb-like, as can be seen from Eq. ( |3.98| ) below: the Coulomb 1/r term is multiplied by 
ln(l/rT), generating a potential singularity in the T — )■ limit. This singularity might be related 
to the instability of our non-BPS string configuration. We discuss how the quantum effects, such 
as string fluctuations, may come in to regulate this logarithmic divergence at very small-r T. 

Using the developed machinery, we also calculate the quark-quark QQ potential in Sec. ^ for 
which there is no hanging (U-shaped) string configuration, and the only contribution comes from 
the exchanges shown in Fig. |14| between the two strings oriented in the same direction. The only 
difference between the QQ potential and the adjoint contribution to the QQ potential is in the sign 
of the NS 2-form contribution, which in the QQ case becomes repulsive. The net QQ potential 
is still attractive, and is plotted in momentum space in Fig. ^ (upper line). It is interesting to 
note that at short distances the QQ potential becomes Coulomb-like, as shown in Eq. ( [4.161 ). The 
absence of the logarithmic singularity that we found earlier in the QQ system can be related to the 
fact that the system of two parallel branes is (BPS) stable. 

We summarize in Sec. |^ by suggesting that the 2-form exchange, which changes sign in going 
from QQ to QQ, may correspond to the chromo-electric modes in the gauge theory, while the 
exchanges of all other supergravity fields (the scalars and the graviton) may be interpreted as being 
mainly due to chromo-magnetic modes. Since both the QQ and QQ potentials are attractive we 
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infer that at strong coupling magnetic modes are more important. We conclude by discussing the 
possibilities of improving on our calculation, including possibly making it more QCD-like. 



2. Two definitions of the heavy quark potential at finite temperature 

Here we will present possible interpretations of the observable we would like to calculate. 
2.1 Singlet and octet potentials 

Imagine immersing a very heavy quark and anti-quark into a static thermal medium of an SU (Nc) 
gauge theory. (For a nice pedagogical presentation of the topic, along with some weak-coupling 
calculations, we refer the reader to [18].) Working in the Euclidean time one defines the Polyakov 
loop operator at a spatial location r by 

L{r) ^Pew(^tg J drAo{f,r)j (2.1) 

with T the Euclidean time, the temporal component of the gauge field, /3 = 1/T, and T the 
temperature. Using Fierz identity one may write 

Tr L\0) Tr L(r) = 2 Tr [i" L^{0) t" L(r)] + ^ Tr [L\Q) L(r)] , (2.2) 

where are the generators of SU (Nc) in the fundamental representation, and f and are spatial 
locations of the quark and the anti-quark. One then defines the color-singlet heavy quark potential 
by 

e-^^^C^) ^ ^(Tr[Lt(0)L(f)]>^, (2.3) 



and the color-adjoint potential by 



e-^^-^W ^ (IV [t^L\0)t^m])^. (2.4) 

Here the angle brackets denote the expectation value of the operators in the thermal bath and the 
matrix elements are normalized to one for the case of no interaction. The distance between the 
quark and the anti-quark is r = |r|. The free energy of the quark-anti-quark pair in a thermal bath 
is given by 

e-'^^'^ = -^(TVLt(0)TVL(r)>^. (2.5) 

c 



2 
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The subscript c in the above equations imphes that we only keep the connected part of the correlator. 
This means that the contributions to the quark-anti-quark free energy and to the potentials due to 
self- interactions of the quarks are subtracted out. (The trivial disconnected part of the correlator 
due to the contribution without any (self-) interactions that gives 1 in the perturbative expansion 
of the correlator is included on both sides of ( |2.5| ).) 

Using the definitions (O), {T§ and (O) in Eq. (P^) we write [17, 18] 



1 



We see that the quark-ant i- quark free energy consists of the singlet and adjoint contributions. One 
can show [18] that at weak coupling and in the large- iVc limit the singlet potential Vi(r) is of the 
order g'^ Nc with g the gauge coupling, that is 



~ A, (2.7) 

A<l,Arc>l 



where X = g^ Nc is the 't Hooft coupling. The singlet potential is also attractive at weak coupling 
[18], and maps onto the standard vacuum Coulomb potential at T = 0. The adjoint (octet) potential 
Vadj{r) is repulsive at weak coupling [18], and is also subleading in the large- iVc hmit, such that 



adj \ 



In the non-perturbative case, Eq. ( p.3|) is used in lattice QCD simulations to determine the singlet 
heavy quark potential at finite-T (see e.g. [26]). 

While the separation of the quark-anti-quark free energy into singlet and adjoint components 
appears to be possible in perturbative calculations order- by-order in the coupling [27,28], it is 
less clear how to accomplish this decomposition in the general non-perturbative case. Indeed the 
definitions of Vi(r) and Vadj{r) in Eqs. ( p.3[ ) and ( p^ are not gauge-invariant. (The sum of the 
two contributions, and, therefore, the free energy, is gauge- invariant, as follows from Eq. (|2.5| ).) We 
hence have a question of whether the singlet and adjoint potential are properly defined. 

One possible resolution of this problem is to define the color-singlet potential using the real 
(Minkowski) time formulation of thermal field theory. Define the temporal Wilson loop W, which 
is a rectangle in the time-r plane, with the spatial sides connecting quark and anti-quark, and the 
temporal sides having length T, as depicted in Fig. |^. The singlet heavy quark potential can then 
be defined by^ 

Vi{r) = lim ^ \nW. (2.9) 

T->-oo 7 



^Clearly this definition of tlie singlet potential is gauge-invariant, and is, in general, different from the definition 



in Eq. (|2.3|). However, the definition (2.9) of Vi{r) does agree with that of Eq. (2.3) both in the lowest-order 



weak-coupling limit and in the strong-coupling AdS/CFT calculations. 
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After this definition of tlie singlet potential, and, after defining the quark-anti-quark free energy 
(in Euclidean time) with the help of Eq. ( p.5|) , one may use Eq. ( p.6|) as the definition of the adjoint 
potential Vadj{r). Below, when calculating the adjoint potential, we will refer to the object defined 
by Eq. Q with the help of Eqs. and (^). 

If one wants to calculate the Polyakov loop correlator ( |2.5| ) (with L 
properly re-defined, as shown in Eq. ( |3.1D below) in the strongly-coupled 
A/" = 4 SYM theory using AdS/CFT correspondence, two distinct string 
configurations have to be considered, as shown in Fig. ^. The first 
AdS/CFT calculations of the finite-T heavy quark potential, performed 
in [8,9] following the zero-temperature calculation of [1,2], found the 
contribution of the hanging string configuration from the left panel of 
Fig. 0. The leading-order contribution of the right panel in Fig. |^ is 
trivial, since, after subtracting the self-interactions of the quarks, which 
correspond to the actions of the two straight strings, we are left with 
zero answer for the renormalized action of this configuration, which is 
what one should expect in the non-interacting case. As we pointed out 
above there are of the order of configurations of straight strings. 
The setup appears to be similar to the decomposition of Eq. ( pIBI) : just 

like in the gauge theory where one has one color-singlet configuration and — 1 color-adjoint 
configurations, there is one hanging string configuration and order- A*"^ straight string configurations 
in Fig. 0. Moreover, the AdS/CFT result for the Polyakov loops correlator can be written as [16] 



n 



Figure 1: Spatial Wilson 
loop described in the text. 



1 



(TrL^(O) TrL(f))^ oc 



-Sha 



9 Jf- ( _ Q-Sstraight 



(2.10) 



with Shanging and Sstraight the actious of the hanging and straight string configurations, with the 
self-interactions removed by renormalization. (Indeed we only know that the coefficient in front of 
the second term on the right of Eq. (|2.1CI| ) is of the order of and we do not have control over 
— 1, which we put there only to normalize the right-hand-side to one in the case of no interactions.) 
Comparing Eqs. ( p.lO|) and ( p. 61) allowed the authors of [16] to suggest that the action of the hanging 
string configuration gives the singlet potential, while the action of the straight strings configuration 
gives the adjoint potential. 

Since, as we have already noted, the decomposition of Eq. ( p.6| ) is not gauge invariant, the 
suggestion of [16] still needs to be proven. However, let us use our real-time definition of the singlet 
potential from Eq. ( ^.91 ). For a Wilson loop from Fig. |1| with a very long temporal extent T, 
the AdS/CFT calculation, carried out in the Lorentzian-signature metric, would only contain the 
hanging string configuration. This is clear since the straight-strings configuration is impossible in 
Lorentzian-signature AdSs: the orientation of the string should be the same throughout the single 
string world-sheet of this configuration, while orientations of the strings connecting to the quark and 
anti-quark are opposite, as shown in the right panel of Fig. |], making this configuration impossible 
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hanging string straight strings 

Figure 2: Two configurations of open strings contributing to Polyakov loop correlator in an AdS/CFT 
calculation. Solid horizontal line denotes the UV boundary of the AdS space, while the dashed line denotes 
the location of the black hole horizon. The arrows on the strings denote their orientations. 



in real time. (In addition to that, in case of Lorentzian-signature black hole the straight-strings 
world-sheet would not be simply connected, indicating additional suppression.) Since the 

Nambu-Goto action of the static hanging string is independent of whether we work in Lorentzian 
or Euclidean signature metrics, the potential V^i(r) defined by Eq. ( p.9| ) is going to be the same as 
found in [8,9]. Thus, using our definition of the singlet potential ( p.9| ), we see that the hanging 
string configuration does indeed give a singlet potential identical to that in Eq. ( p.3| ), yielding [8,9] 



Vl(r) 



~ \/A. (2.11) 

Afc>A>l 



The fact that this potential is of the same order in Nc as the weakly-coupled singlet potential in 
Eq. ( p.7| ) indicates consistency of our conclusions with perturbative calculations. 

Using Eq. ( |2.6| ), and remembering again that the hanging string configuration gives the singlet 
potential regardless of whether we work in real or imaginary time, we see that the straight strings 
configuration from Fig. ^ gives the adjoint potential. At the leading- A*"^ order without any interac- 
tions between the strings the adjoint potential defined this way is zero. As was suggested in [11], 
a non-zero contribution to Sstraight comes from the interactions of the strings due to exchanges of 
supergravity fields, as shown in Fig. ^. Since the coupling of the supergravity fields to the string 
world-sheet is of the order of 1/A"c, the contribution to the two-strings action due to the exchanges 
shown in Fig. |^ is order-l/A^^. Exchanges of supergravity fields between various string configura- 
tions have been considered before [29-31], yielding the interaction of the order of X/N^. Identifying 
the diagram in Fig. ^ as the leading non-trivial contribution to the adjoint heavy quark potential, 
we conclude that at strong coupling 



adj \ 



- W^- (2.12) 
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Comparing this result to Eq. ( p.8|) we observe that, surprisingly, for the adjoint potential the power 
of 't Hooft coupling A does not change in going from small to large A! This is unlike the case of 
the singlet potential in which, just like in many other AdS/CFT results, the power of the coupling 
changes from A to as we increase it from very small to very large, as can be seen from Eqs. ( |2.7| ) 
and pni ). 

In the terminology of the two heavy quark potentials defined 
above, the goal of this paper is to calculate the contribution to the 
adjoint potential depicted in Fig. |. 

2.2 One unified potential 

There exists an alternative to the singlet and adjoint potential 
decomposition.^ One can simply define one unified heavy quark 
potential V{r) as the free energy of the quark-anti-quark pair in 
the plasma, such that [32] 

1 

m 




-l3V{r) 



(TrL^(O) Tr L(f))^ 



Figure 3: A correction to 
(2.13) the straight-strings configuration 
from Fig. ^ due to exchange of a 
The strength of this definition is that the heavy quark potential supergravity field, which, as we 
defined this way is manifestly gauge- invariant. 

However, it appears difficult to find an intuitive physical in- 
terpretation for the single unified potential defined in Eq. (|2.13|) . 
To see this, let us start with Eqs. (pD and (P), which, together ^^o^^^ "^^^^ a cork-screw hne. 
with Eq. ( p.l3| ) give 

1 



will see below, can be either a 
dilaton, massive scalar, 2-form, 
or a graviton, all of which are de- 



-/3V{r) 



-l3Vi{r) 



]■ 



(2.14) 



Remembering from the previous Section that Vi(r) ~ A*"^ and Vadj{r) ~ l/A"^ (both at small and 
large coupling), we expand the right- hand- side of Eq. ( p.l4| ) in powers of to obtain [11] 



-fiVir) 



-/9Vi(r) 



1] 



We see that the unified potential V{r) has to be iV^-suppressed and equal to 

1 1 



V{r) 



[1 



-/3Vi(r) 



(2.15) 



(2.16) 



It seems a bit counter-intuitive that the heavy quark potential at T 7^ should be A^^-suppressed, 
while the T = potential is not. Also, in the T — )• limit, the singlet potential Vi{r) maps 
smoothly onto the vacuum potential both at small [18] and large [1,2] couplings: it is clear from 



^We would like to thank Larry Yaffe for a discussion on this subject. 
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Eq. ( 2.16 ) that the potential V{r) instead goes to infinity in the T — )■ hmit. Moreover, at small 
quark-anti-quark separations r, both in QCD and in N' = 4 SYM theory, the singlet potential is 

Vi{r) 1/r (while Vadj{r) ~ 1/r in QCD [18], and Vadj{r) (1/r) ln(l/r T) in AT = 4 SYM 

at large- A, as we will show below). However, as one can infer from Eq. ( p.l6| ), the potential V{r) 
diverges exponentially at small-r as — exp(const /3/r), i.e., it becomes very strongly attractive: the 
physical origin of this behavior is not clear. 

Despite the problems with its interpretation, the potential V{r) is well-defined. As we have 
already mentioned, in our calculation below we will find the contribution of the exchanges of su- 
pergravity fields between the two straight strings pictured in Fig. |^. In the language of the unified 
potential V{r) in Eq. (p.l6|) , we will be constructing one of the contributions to this potential 
(Vadjir)) at order-l/A"^, while the other contribution (coming from Vi(r)) was previously found 
in [8,9]. 



3. The adjoint contribution to the QQ potential 

We want to find the contribution to the correlator of two Polyakov loops (TrL^(O) Tr L(r^)^ in 
A/" = 4 SYM theory at strong 't Hooft coupling at order- l/A*"^. On the string theory or gravity 
side, these contribution will be coming from the interactions between the two straight strings in 
Fig. ^ . In Af = 4 SYM theory the definition of the Polyakov loop operator suitable for AdS/CFT 
calculations is different from Eq. ( |2.1| ), and is given by [1,29] (in Euclidean time) 

L^=4(rl = Pexp (^g j dr [Ao{f,T)-i9\f,T)X\r,T)]^ , (3.1) 

where X\f, r) with I = 1, . . . , 6, are the six scalar fields of A/" = 4 SYM, and 6^ is a point on the 
unit five-sphere, which, for simplicity, will be taken to be independent of time. To make our results 
more QCD-like we will try to reduce their dependence on the scalar fields by averaging Lj^=4 
over all angles 9^ . 

At low energies type IIB string theory compactified on AdSr, x contains a spectrum of 
supergravity modes. In the bosonic sector, the lightest modes are graviton Hmn, dilaton 0, axion 
Co, Neveu-Schwarz (NS) Bmn and Ramond-Ramond (RR) Cmn 2-form fields (M, A^ = 1, . . . , 5), 
and RR 4-form field C4 with a self-dual field strength. The spectrum also contains an infinite 
tower of massive scalars, t'^ {k > 0) and {k > 2) with k the index of KK modes [33-35]. We 
are interested in those modes that couple to a string world-sheet (with the coupling that is least 
suppressed in Nc). The smearing of the string over the will effectively eliminate the coupling of 
the string world-sheet to modes that carry non-zero KK charge. In what follows, we will neglect 
the contribution of the fermionic sector to the interaction of the two string world-sheets, since it is 
A'^-suppressed compared to that of the bosonic fields. 
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We want to find the quadratic fiuctuations of those fields that couple to the string world-sheet 
as shown in Fig. |^ in the background of the AdS Schwarzschild black hole (AdSSBH) metric 



f{z) dr"^ + dx'dx' 



dz' 



1-^ 



z 



(3.2) 



h 



where z G {0, Zh) is the direction along the extra fifth dimension, = (x ,x , x ) are spacial 
coordinates, r is the Euchdean time with period (3 = Tczh, and z^ determines the position of the 
black hole horizon, that is related to temperature as follows: T = = l/lnzh). 

To find the correlator of Polyakov loops at order- we will only need to consider exchanges 
of the fields that couple to string world-sheet at the tree level. Such fields are graviton, dilaton, 
massive scalars, and NS/RR two forms. The RR 2- form does not couple directly to a string world- 
sheet but it mixes with NS-NS 2-form, making it effectively massive [15,33]. Notice that at the 
tree-level the axion only couples to the world-sheet fermions, which leads to a higher suppression 
in Nc- Since we are interested in order-l/A^^ corrections, we will ignore the contribution to the 
correlator coming from the exchange of the axion. 

As we argued above, the quantity we would like to calculate is 



TrLt(0,^') Tr L(f,^) 



TrLt(0,^')) (TrL(f 



straight 



adjoint 



S straight ~^ ^ ( ^4 



(3.3) 



where we assume that the quark and anti-quark Polyakov line operators ( |3.1| ) are taken at fixed 
angles 6 and 6' correspondingly. Using the results of [29] we can write this contribution (somewhat 
schematically) in terms of the integrals over string world-sheets A and A' 



TrLt(0,^') Tr L(f,^) 



TrLt(0,^')) (TrL(f,^) 



exp 



adjoint 



^ — ' / Zn a Zn a 

. k,m 



(3.4) 



Here Gk{z,T; z' ,t') is the sum of bulk-to-bulk propagators (multiplied by the appropriate vertex 
factors different for each field) for all the supergravity fields that couple to string world-sheets 
(parametrized by z,t and z',t'), with the appropriate indices chosen for the graviton and 2-form 



contributions. The functions Ym '^ (6) are spherical harmonics, Ym' (6) = Cf^ Oi^ - ■ ■ ^4 with 
labeling all harmonics corresponding to the state with the total angular momentum = k{k + 4) 
and C^l a basis of symmetric traceless tensors such that C]^^ (jh---ik ^2 _ §mim2 [29^36]. As 
usual the slope parameter a' = 1/ \/X. 

Expanding Eq. (^) to the first non-trivial order, and averaging it along with Eq. ( p.3| ) over 
angles 6 and 6' yields 



m 



^ {S straight) 



1 



dA dA! 



2 vr a' 2 vr a' 



Gniz, r; z\ t') 



(3.5) 
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since only k = KK mode survives the averaging and Yq = 1 in the conventions of [29,36] that we 
have adopted. Below we will calculate k = KK mode contributions to the action Sstraight coming 
from the dilaton, massive scalar 2- form fields, and the graviton. 

3.1 The Dilaton 

3.1.1 Dilaton potential and the EOM 

We begin with the simplest case of the dilaton exchange between the sting world-sheets pictured 
in Fig. ^. Writing the KK expansion for the dilaton field as </> = (pkY^^^O), and only taking the 
lowest harmonic, the 10-dimensional dilaton action will be reduced to 



Ar2 

Q C 



16 vr^ 



j d'^x^ g^'''dM<pdN<P , (3.6) 



where = 0(x'^, z) is the coefficient of the zeroth KK harmonic and M, N = fi,z with /i = 0, . . . , 3. 
Notice, that the massless dilaton in the AdS bulk is dual to the glueball operator TrF^^ in the 
boundary theory. The coupling of the dilaton to the string world-sheet is described by the action 

^— J rfVe^/2V7 + ^ / d'a^R^'U, (3.7) 

where ja/s = QhiNdaX^ dpx^ is the induced metric on the string world-sheet, and i?^^^ is the world- 
sheet curvature. In the large N,. limit, when the AdS radius is large, R^"^^ will be subleading and 
can be ignored [29]. 

Expanding e'^^'^ = 1 + 0/2 + . . . and using a' = 1/VX we see that the leading-order action of 
the dilaton coupled to two strings is 

S^^ = j d''x^g''^dM<pdN<P + ^ I d'a^(l) + ^ [ d'a^cj), (3.8) 

16 TT^ J An 7(1) An J ^2) 

where /(^^ has a meaning of integration over the i^^ string world-sheet, i = 1,2. The equations of 
motion (EOM) for the rescaled scalar field, 

^ ' (3.9) 



27r 



can be written as 



[^9 9''''^N^ = ^5^'\x-X,) , (3.10) 

where sub- index i = 1,2 implies that we want to find the scalar field created by the string localized 
at spatial position Xj. Choosing the quark to be located at Xi = -Fand the anti-quark at X2 = 0, 
the two strings world-sheets can be parametrized as 

Xf = (r, f, z = a), Xf = (r, 6,z = Zh-a) (3.11) 
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with r G [0,/3] and a G [0, z/j]. Note that in Eq. the string dual to the anti-quark (X2) is 

oriented opposite to that of the string dual to the quark (Xi), as also shown in Fig. ^ this will be 
important for the analysis of the 2-form fields contribution. 

Now, we are interested in evaluating the action ( |3.8|) on the solution of the classical EOM 
( p.lOp . Taking into account that = 0i + 02, and dropping the terms in the action that correspond 
to (anti-) quark self-energy contributions as they do not contribute to the connected Polyakov loop 
correlator in Eq. (|3.3|) that we want to find, we obtain 



S, 



dil 



i TT 



— / d 02 + ^ / d (y^/l 



(1) 



^1 , 



(3.12) 



(2) 



where integration in each term is over one of the string world-sheets of the classical field created 
by the other string. Note that the two terms on the right-hand-side of Eq. ( |3.12| ) are equal. The 
integration over string world-sheet is 



dr J dz. 




(3.13) 



The determinant of the induced metric is 7 = det'^ap = for the strings parametrized as in 
Eq. ( p. lip in the background metric ( p.2| ). Using all this in Eq. ( |3.12| ), and noticing that in the 
static case considered here the solution of classical EOM ( |3.10| ) is going to be time-independent, we 
obtain the dilaton contribution to the potential 



A 



2 m 



(3.14) 



As 02 is time- independent, we explicitly removed r from its argument. 
Since g = det qmn = 1/^^° the EOM (|310D for 02 is 



f 



The equation is easier to tackle in momentum space. Writing 

d^q 



h[r,z) 



{2n) 



■e^^Xg,^) 



we recast Eq. ( p.l5|) as 



z^\ 1 



1-^ ^^^- 3 + ^ -fz-q f = z 



4/ z 



(3.15) 



(3.16) 



(3.17) 
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with = q"^, ip^^ = dl^p, ip^ = d^fp. In solving Eq. ( |3.17|) it is convenient to measure z and q in 
units of Zh, such that, rescahng 

> z, qzh q, ^ (3.18) 

to make them dimensionless, we get 

(1 - z^) (^,, -[^ + z^)-^,-q^ip = z. (3.19) 

z 

The potential in momentum space Vadj (q) is defined by 

K.,(r) = - / 7^e^^^-^>ad,(g), (3.20) 
Zh J [2 Try 

where we have switched to dimensionless f/zh — ?■ r as well. The dimensionless momentum-space 
dilaton contribution to the potential is 

1 

KV^) = ^ / $^(^,-)- (3.21) 



(Note that the rescaled (p in Eq. ( p.21| ) is dimensionless, and so are the rescaled z and q.) Since 
Eq. ( |3.19| ) depends on q through g^, its solution ip depends on the vector q through its length 
squared q^. 

3.1.2 Solution of the dilaton EOM 



Solution of the dilaton EOM ( p.l9[ ) can be constructed as a series in the powers of z 



oo 



h ^2n+3 °° f ^2n+4 

n=0 n=0 



with the recursion relations for the coefficients 



q^ (f' 

bn = -TZ TTTT^ + ^"-2, &0 = -1, h = - — , (3.23a) 

[2 n — Ij [2 n + 1) 6 



fn — 777 —7^7 fn-l + fn-2, fo — 1; /l — (3.23b) 

2n{2n + 2j 8 

The first series in Eq. ( p.22| ) is the solution of the inhomogeneous equation, while the second 
series is the solution of the homogeneous equation. Both series were constructed by requiring that 
(p{z = 0) = 0: this is the condition of normalizability of the dilaton field. Clearly the potential 
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in ( p.21 ) would be infinite if tliis condition was not satisfied. Tlie coefficient C'^(g^) in Eq. ( p.22 ) 



has to be fixed by the boundary condition at the horizon. To obtain this condition we follow the 
standard procedure and expand the AdSSBH metric ( |3.2| ) near horizon. Defining p = 2y/l — z we 
obtain a fiat near-horizon metric 

In this new flat metric the function (p{x, p) should have zero derivative at the origin, <y3p(p = 0) = 0, 
otherwise ip would have a discontinuous derivative at the origin.^ This translates into 

{Vl^z^,} 1^^^ = (3.25) 



in z-space. The condition ( p.25| ) appears to be somewhat weak, as it only implies that ip should be 



finite at the horizon. To see that the finiteness of if at the horizon uniquely determines the solution 
of Eq. ( |3.19| ) we can expand (p in the powers of 1 — z. In this near-horizon expansion case the series 
solving the inhomogeneous and homogeneous equations ( |3.19| ) are the same, such that the solution 
can be written as a single series 

oo 

if = Y,hn{l- zY (3.26) 

n=0 

with the recursion relations 

_2(n-l)(5n-8) + g\ 2 (n - 2) (5 n - 12) + 

hn — -. 5 hn-l 5 hn-2 

4 4 



(n-3)(5n-16) ^ (n - 4)\ 1 ^^A 97^ 

4:71'^ 4 8 3d 

between its coefficients. One can see that the value of at the horizon z = 1, given by the coefficient 
ho, is assumed to be finite here. The coefficient ho alone determines the rest of the series (|3.26|) , 
with the help of Eq. (|3.27|) . The only remaining degree of freedom, ho, is fixed by the ip{z = 0) = 



condition. This proves that finiteness of (f at the horizon, along with the (p{z = 0) = condition, 
are sufficient to uniquely define the solution of Eq. ( p.l9| ) . 

Analytic summation of the series ( |3.22| ) and/or ( ^.2(j| ) appears to be prohibitively complicated. 
Instead we construct the solution for ip by summing the series ( |3.22D numerically. We numerically 
evaluate partial sums 

^{q, z,N) = J2 + ^^^^ N) J2 (3-28) 

n=0 n=0 

constructing sequential approximations ip(q,z,N) of the exact solution ip(q,z) = lim ip(q,z,N). 
To determine the coefficients C'^{q'^, N) we need to require that ip{q, z, N) is finite at the horizon: to 



^Wc thank Samir Mathur for pointing out this argument to us. 
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insure this without specifying the (unknown) value of ip{q, N) at the horizon we impose Neumann 
boundary conditions on ip{q, N) aX z = 1. Indeed the exact solution of Eq. (p.l9|) does not have to 
satisfy Neumann boundary conditions at the horizon: however, imposing this condition on partial 
sums ip{q, z, A^) makes them finite at the horizon, such that, since the solution which is finite at the 
horizon is unique, as increases the partial sums would converge onto this one exact solution for 
z), with the part of the interval z G (0, 1) affected by the Neumann boundary condition rapidly 
shrinking. Since Eq. ( p.l9| ) is a second-order differential equation, it has two solutions: one that 
is finite at the horizon, and one that is infinite at the horizon. Therefore any boundary condition 
forcing (p{q,z, N) to be finite at 2; = 1 would map these partial sums onto the solution finite at the 
horizon in the N ^ 00 limit. For instance, Dirichlet boundary condition at z = 1 would also work. 
However, Neumann boundary conditions appear to give the fastest numerical convergence to the 
exact solution. 



'P(l,z,N) 



-0.02 



-0.04 - 



-0.06 



-0.08 - 



-0.10 




Figure 4: Numerical evaluations of (p{q, z, N) as functions of z (in units of Zh) from Eq. (|3.2^ ) for q = 1 
and = 1, 3, 5, 11, 15, 21, 31 in the descending curve order. 



Demanding that dzip{q, z, N) 



z=l 



yields 



C^q',N) 



Yln=0 



(3.29) 



Using Eq. ( p.29|) in Eq. (|3.28|) we construct a series of numerical approximations to the solution of 
Eq. ( p.l9| ). Three approximate solutions are plotted in Fig. ^ for different values of A^. One can 
see that, as the order of the partial sum A^ increases, the solutions tend to converge to the unified 
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limiting curve, deviating from it only in the rapidly shrinking region near z = 1 to satisfy Neumann 
boundary condition. 




Re{q) 



-0.10 



Figure 5: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the dilaton field between the string world-sheets, plotted along the real-g axis (Im g = 0) in units of 
\/{2N^). The plot is for N = 10^ iterations. (Here and in all subsequent plots we put Zh = 1.) 

Since we are interested in the dilaton contribution to the heavy quark potential V^^j{q) we use 
Eq. { ^2^ in Eq. {^2^ to obtain 



A 



2 AT 2 



n=0 



(2n + 2) (2n + 3) 



fn 



oo 

2 (2n + 3) (2n + 4) 



(3.30) 



The potential is also evaluated by numerical calculation of the partial sums. Dilaton contribution 
Kidj (q) ^'^ heavy quark potential is plotted in momentum space in Fig. ^ The potential is finite 
at g = which indicates screening of the quark-anti-quark interactions by the thermal J\f = 4 SYM 
medium. Since the potential is monotonically increasing and negative, one should expect that in 
coordinate space the contribution of the dilaton will be attractive. 

3.1.3 Asymptotics of the dilaton contribution 

While it is hard to perform a reliable numerical Fourier transform of the potential from Fig. |^ into 
coordinate space, we can use our numerical results to understand its asymptotics at large and small 
values of the dimensionless parameter r, or, equivalently, of the parameter rT if we go back to the 
dimensionful r not measured in the units of Zh- 
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Let us begin with the case of large-r T, which can be achieved by either increasing the quark- 
anti-quark separation r or by increasing the temperature T. Noting that the momentum space 
potential is a function of we can integrate Eq. (|3.20|) over the angles, obtaining 



oo 

K.,(r) = J^,^^ I dqqe^^^V^,M')- (3-31) 

— oo 

The g-integral can be done by closing the contour in the upper half-plane and picking up contri- 
butions of all the singularities there. It turns out that the singularities of the dilaton potential 
V^^-{(l) (along with the similar potentials for other supergravity fields that we will analyze below) 
are only limited to poles along the imaginary-g axis. As originally suggested in [11], the positions 
of the poles correspond to the glueball masses for QCD3 calculated in [14, 15,37-40] and references 
therein. The potential V^^j{q) is plotted in units of \/{2N^) along the positive imaginary-g axis in 
Fig. ^ illustrating this point. Note that the potential is real along both the real and the imaginary 
axes in the complex g-plane. 

Denoting the positions of the poles on the positive imaginary-g axis by zm„ with n = 1,2, . . . 
(such that mi < m2 < . . .) we write (see also [30,37])'' 

^ 00 

Vadjir) = V^m^e-'""'- lim [(g - ^ m„) Kd,(g')] • (3.32) 

" n=l 

The large-r behavior is determined by the first pole mi. 

For the case of the dilaton, according to our numerical solution the leading pole is at mf = 
3.4041 ± 0.0001 (for = 3 x 10^ iterations), in agreement with [14, 15]. The residue of the pole is 
i (0.12 ± 0.01), such that the large-r asymptotics of the dilaton contribution to the potential is 



A 0.41 ±0.04 ^.3,0.... (3 33) 



where we have inserted Zh = 1/{tt T) back into the expression. We obtain a screened Yukawa- type 
attractive potential falling off exponentially with the distance r. 

Now we consider the case of small-r, which, inserting Zh back, means small rT, a regime which 
can be interpreted as resulting from either the short separations r or low temperatures T. Due to 
the Fourier transform ( p.20| ), low-r corresponds to large-g. In the large-g limit, and for z ^ 1/g, 
we can neglect the terms on the left-hand- side of Eq. (|3.19|) containing 2;-derivatives compared to 
the term with g^. This yields 



^In Eq. (3.32) we assume that the poles are order one: we have exphcitly verified this only for the first pole for 



each supergravity particle analyzed. 
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Figure 6: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the dilaton field between the string world-sheets, plotted along the imaginary-g axis (Re g = 0) in units 
of X/{2N^). The plot is for N = 5000 iterations. 



Since, as follows from the series solution ( ^.221 ), the function if goes to zero as when z — )• 0, we 



can argue that most of the support of the integrand (p/z"^ in Eq. ( p.21 ) comes from larger z, and 
hence, for large-g, the condition g^; ^ 1 is satisfied over most of the z-range contributing to the 
potential. Using Eq. ( |3.34|) in Eq. ( p.21|) and inserting a UV cutoff of 1/g in it, we obtain (with the 
leading logarithmic accuracy) 



2iV2 q 



(3.35) 



Fourier-transforming this back into coordinate space with the help of Eq. (|3.20| ), and inserting z^ 
back in to make r dimensionful yields (again with the leading logarithmic accuracy) 



X 



rT<l 



— In — 

2N? 47rr ^ rT' 



(3.36) 



The potential in ( p. 36] ) is attractive, but is not quite Coulomb-like, due to the extra logarithmic 
factor. This factor is somewhat worrisome, as it makes the potential infinite in the T — )■ limit. 
This, however, might not be viewed as a problem on the gauge theory side, since it is not clear 
how to define the adjoint potential at T = in a gauge- invariant way: therefore, the divergence 
takes place in the limit when the quantity we are calculating is not defined. On the string theory 
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side it appears that, in absence of the black hole, the interaction between two string world-sheets 
would be infinite: the situation is analogous to the instability of the parallel D-brane-anti-D-brane 
configuration. One may also argue that at small separations r higher-order string corrections may 
come in and modify the result, possibly regulating the divergence: we will discuss this possibility 
after calculating the contributions of all other supergravity fields. 

3.2 Massive Scalars 

3.2.1 The potential and EOM 

In addition to the dilaton, compactification on 5*^ leads to additional scalar modes that come from 
the KK modes of the trace of the metric and from the 4-form field perturbations over S^, that is, 
from /i" and ttap^s correspondingly, where indices take values on [29,33-36]. The KK modes 
of h^, and aap-ys mix in the EOM, but at the quadratic order the action can be diagonalized by 
the scalar modes s'^ (with k > 2) and t'^ (with k > 0) [29,33-36]. The modes for k = 2,3 are 
tachyonic, while all other s'' and all t'^ modes are massive: we will refer to all modes t'' and as 
massive scalars. Since we are interested in the modes which survive averaging over S^, i.e., in k = 
KK modes, we only need to study the contribution of to the string interaction. The action for 
the modes can be found in [34,35]. The coupling of to the string world-sheet follows from the 
coupling of the graviton to the string. Writing the 5d metric fiuctuations for k = KK mode in 
terms of fields diagonalizing EOM as 



^0 

^9mN = ^MN ^dMNt + -^'^Al'^Nt (3.37) 

with Vm the covariant derivative we can determine the coupling of t° to the string world-sheet from 
the couphng of Sg^jj^ to the string. One can readily show that the coupling of the Vm VArt° term 
to the world-sheet of our straight strings (|3.11|) is zero if dzt^{z = 0) = 0, with the latter condition 
satisfied by the solution of the classical EOM for Dropping the derivative term in Eq. ( p. 371 ) we 
get (using the same notation as in Eq. 



o 7X J o 7t J Z o 7C J Z 

(1) (2) 

(3.38 



Rescaling the t° field 



p^Si^^o (3 39) 



we write the equation of motion for the field produced by the string at X2 = 0: 



^ Om [v^^7^^^ d^t 0] - 32 1 ° = (x) . (3.40) 

V 9 
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The action evaluated at the classical solution is 



c 





6 N"^ 

giving the contribution to the heavy quark potential 



/3 

St'> = l4-Jdr f^^t^r^z) (3.41) 



Vfjr) = 14 / ^i\r,z). (3.42) 





Just like with the dilaton, in order to solve Eq. ( p.40|) we go to momentum space 



(2vr) 



where, in units of Zh ( p.l8| ), the EOM becomes 



(l-.^)tL-(3 + .^)lt°-g^t°-§t° = .. (3.44) 



z ' z'^ 



Just like with the dilaton we require that t°(2; = 0) = and t°(2; = 1) is finite. The contribution of 
t° to the momentum-space potential is then (in the same Zh = \ units) 



1 



V» = ~/$l°(«.^). (3.45) 

^ 



The coordinate-space potential can be obtained from Eq. (|3.20|) . 

3.2.2 Solution of the EOM for t° 

The series solution for Eq. ( |3.44| ) is 



oo 



n=0 71=0 

with the recursion relations for the coefficients 

_ {2n-\f _ 1 _ g2 

" (2n-5) (2n + 7) """^ +(2n-5) (2n + 7)""-'' "° " ~35' " 945' ^^'^^^^ 

, (n + 2)^ g2 
4 n (n -I- 6) n (n + 6) 28 



- 20 - 



2N}.Vi 




-0.10 



Re(?) 



Figure 7: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the field between the string world-sheets, plotted along the real-g axis (Im g = 0) in units of A/ (2 N'^). 
The plot is for = 10^ iterations. 



Again we sum the series ( p.46|) numerically, imposing the Neumann boundary conditions at the 
horizon. These yield 



En=0 Cn{2n + 3) 



(3.48) 



Substituting Eq. ( |3.46| ) into Eq. ( |3.45| ) we obtain the contribution of the scalar t° to the 
momentum-space heavy quark potential 



6 m 



n=0 



9n 



n=0 



2n + 7 



(3.49) 



The potential V^lj{q) is plotted in Fig. in units of X/{2N^). It looks similar to the dilaton 
contribution to the potential plotted in Fig. |^ and also appears to be attractive. 

3.2.3 Asymptotics of the massive scalar contribution 

Just like for the dilaton, the singularities of the potential V^dji^) poles along the imaginary-g 
axis. The potential V^'^j (g) along the positive imaginary-g axis is plotted in Fig. ||. 

The large-r asymptotics is given by the leading pole, which was found to be at m\ = 7.410 ± 
0.001 (for = 10^ iterations), in minor disagreement with [15], possibly due to differences in 
implementation of the numerical simulations. The residue of the pole is i (0.04 ± 0.01), giving the 
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Figure 8: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the field between the string world-sheets, plotted along the imaginary-g axis (Re g = 0) in units of 
X/{2N^). The plot is for N = 10^ iterations. 

large-r asymptotics of the t° contribution to the heavy quark potential 



with Zh = 1/{ttT) inserted back such that r has units of distance. Again, this is an attractive 
screened Yukawa-type potential, with the screening mass larger than that for the dilaton. Hence 
the dilaton contribution dominates over that of at large-r T. 

We have also investigated the contribution of the tachyonic scalar = s'^^^ to the interaction 
between the string world-sheets. Since this is a A; = 2 KK mode, it does not contribute to the 
potential ( p.5|) averaged over positions of both strings that we want to calculate. In the effort to 
translate the physical meaning of the contribution into the QCD language, and identifying the 
SU{4:) i?-symmetry group of A/" = 4 SYM with the flavor group on QCD, one may argue that s^- 
exchange would correspond to some flavor- changing interaction between the quarks. Our motivation 
for studying lies in the fact it is a tachyon, and it is precisely the contribution which dominates 
the large-r asymptotics of the interactions of two strings in the empty AdSs space (zero-temperature 
case in the gauge theory), as was shown in [29], giving a potential between two rectangular Wilson 
loops (two "mesons") that falls of as 1/r^ at large-r. The s^-exchange in the AdSSBH background 
leads to an attractive potential, but with a small screening mass mf ^ 1.373 ± 0.001 (for = 10^ 




rT>l 



5 A 0.30 ±0.08 

6 Af 2^7 



(3.50) 



- 22 - 



iterations), such that 

A 1 



oc-^^e---'^- (3.51) 

T>1 ^^c ' 



Clearly, if we allow for exchanges of higher KK modes, the contribution would dominate over that 
of the dilaton and at large-r T. We would like to point out in advance that the screening masses 
for the 2- form field and the graviton (A; = term only) would also be larger than mf ~ 1.373. 
As the screening mass appears to increase with fc, one may conjecture that the contribution 
( p.51| ) would dominate over all other supergravity modes exchanges at large-r T. Since our main 
analysis here is limited to = case when quarks do not carry i?-charge, we will not have the 
contribution in the net potential. 

At small-r/large-g the approximate solution of Eq. ( |3.44| ) is 



giving the coordinate-space potential 



'\.». - (3-52) 



where again we have reinstated the temperature T explicitly. 

3.3 Two-Form Fields 

3.3.1 The potential and EOM 

The effective action for the A; = KK mode of the NS 2-form along with its coupling to the string 
world-sheet can be written as (after integrating out the contribution of the RR 2-form) 



24/^2 J --v^"-.v«^^ , ^2 



+ w / ^''^ ^"^"^ ^'^"^ + iJ^ / ^'^^'^ ^'^"^ ' ^^-^^^ 

(1) (2) 

where indices a, & = 0, 1 run over string world-sheets, e"'^ = — e^" with = 1, and = An^/N^. 
As usual 

Hmnr = Qm Bnr + On Brm + Qr Bmn- (3.55) 

The coefficient ^ is equal to 1 when both indices M,N of Bmn in the same term are spatial, i.e., 
when M, N = fi,!/ = 1, . . . ,4, and = i when either M or is 0. This extra factor of i arises due 
to the transition to Euclidean-signature space with ^Mmkowskmn _ ^^Euchdean^ where fx = 1, ... ,4 
now. 
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The second term on the right of Eq. ( 3.54|) describes the effective mass of the 2-form that is 
generated via the Chern-Simons type of mixing with the RR 2-form Cmn [15,33]. To understand 
why there should be a mass term, notice that, in the absence of the string sources, the 2-forms 
Bmn and Cmn can be combined into one complex 2-form field Amn = Bmn + iCMN [15,33]. The 
field equation for Amn can be factorized into two first order differential equations, and each can be 
iterated leading to a second order equation of the form Max Amn — Amn = 0, where Max is a 
second order Maxwell differential operator for 2-forms defined by Max Bmn = Hmnr (for niore 
details see [15]). For = KK modes that have no dependence on coordinates, one equation 
has m? = 16, and the other one has m = 0. It can be shown that the massless equation has only 
pure gauge solutions, which can be ignored (see [33]), leaving us with the massive equation for the 
combined field Amn- Taking real and imaginary parts of this equation (working for the moment 
in the Lorentzian-signature metric where both Bmn and Cmn are real) we would get separate 
differential equation for Bmn and Cmni namely Max Bmn~^Q Bmn = and Max Cmn — Cmn = 
0. Since only the NS 2-form field Bmn couples to string world-sheets, we are interested in the 
resulting action for it that gives this massive EOM, which is given by the first two terms on the 
right of (13:5^) . 

Using the string parametrization (|3.11|) in the action ( |3.54|) simplifies it to 

- J d^'x^ Hmnr H''''^ + ^ y ^'^v^ 



S 



B 



MN 



96 7r2 



+ 



27r 



dr dz Bqz{t, f, z) 



27r 



dr dz B, 



02 



T,0,Z), 



(3.56) 



(1) 



(2) 



with the action along the classical solution 

O R 



2h 



2n 



dr / dzBoJT,r,z) 



(3.57) 



To find the classical field Bmn{j^ z) due to the string X2 we define the rescaled 2-form field 



B 



MN 



B 



MN- 



47r 

We need to find the classical field Bmn satisfying the EOM which follow from Eq. (|3.54| ) 



(3.58) 



—p^ 9mm' 9nn' dji 
V9 



16 5 



MN 



(5'Mo Qnz - Ono Qmz) - (3.59) 



The solution of this equation, along with Eq. ( p. 571) , would give us the contribution of the 2-form 
field to the adjoint QQ potential 



2 \ r 

— J dzBoz{r,z). 



(3.60) 
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To find Bqz we need to solve the EOM (|3.59|) to find the field of the string X2. In the static 
case we consider, and due to rotational 0(3) symmetry around the string in the i?^-space spanned 
by x^,x'^,x^, the solution should depend on r and z only, Bmn = BMN{r,z). Working in spherical 
coordinates r,6,(f) in we notice that EOM for Bmn components with either (or both) M or N 
equal 6 or (p decouple, and since there is no source for those components we can put them all to 
zero. This leaves us with Bq^, Bq^, and Brz- The r z component of Eq. ( |3.59| ) requires that B^z = 0, 
leaving us with Bq^ and i?or only. 

The EOM for components Bqz and i?or are mixed with each other. The z and r components 
of Eq. (|3.59|) yield (keeping in mind that the fields depend only on (r, 2;)) 

(0^) ^dr[r^drB,z-dzBor)]-lQBoz = z5%r} , (3.61a) 
(0 r) z f{z) dz [z {dr Boz - dz 5o,)] + 16 So. = . (3.61b) 
Defining A(r, z) by dr A = i?or we can shift the 2- form field 



Oz 



O2 



and use the new field Bqz to simplify Eqs. ( |3.61| ) to 

r'^drBoz 



zB, 



Oz 



zf{z)d. 

Fourier-transforming the 2-form field 

BQz{r,z) 

and eliminating A from Eqs. ( p.63|) yields 



a A 



-16 [B^z + dzA 
+ 16A = . 



z5'^{r) 



(2^ 



zH\ 



El, + ;z (3 - 7;z^) E'^, - (15 + 5 + z^) E, 



Oz 



(3.62) 

(3.63a) 
(3.63b) 

(3.64) 

(3.65) 



Here the prime denotes the partial derivative with respect to z. We have also switched to the units 
of Zh (see Eq. ( p.l8|) ) in Eq. ( p.65|) with Eoz/zh — )■ Eqz- The boundary conditions for Eqz are the 
same as before: Eoz{z = 0) = and Eoz{z = 1) should be finite. 

Note that Eq. ( ^ISBj) insures that A{z = 0) = A{z = Zh) = 0. Together with Eq. ( ^I^ 
these conditions allow one to replace Bqz in Eq. (|3.6CI|) with Bqz without changing the value of the 
integral. In momentum space the contribution of the 2-form field to the adjoint QQ potential (with 
the Fourier transform defined by Eq. ( |3.20 )) is 

1 

K|(9) = 1^ / dzEoziq,z) (3.66) 



with z and q now taken in the units of Zh- 
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3.3.2 Solution of the 2-form EOM 

The series solution of Eq. ( ^.651 ) is slightly more involved than that for the dilaton and t° and is 
given by 



n=0 n=l n=0 

with the recursion relations 

i2 



2n+3 



dr. 



A{n-l){n + 3) 



"n-l H —^T "n-2 

n + 6 



n + l 



■ln + 



Z 



n 



(3.67) 



{n -l){n + 3) (n - I) {n + 3) 



n-2; 



4(n- 1) (n + 2) 



n + 3 



'n-25 



n + 2 
■ e„_i H ; — 7 e„_2, 



/o = 0, h 



e_i = 0, Co 



96^ 



4n(n + 4) ' ' n + 4 
Neumann boundary conditions for the partial sums of the series give 

dn{2n + l) + Y.n=l In 



T.n=Q e„ (2 n + 3) 



(3.68a) 
(3.68b) 

(3.68c) 
(3.69) 



Plugging Eq. ( p. 671) into Eq. (|3.66|) we obtain a series representation of the contribution of the 
2-form field to the heavy quark potential in momentum space 



2A 



dr. 



2n + 2 



2n + 4 



(3.70) 



_n=0 n=l ^ ' 71=0 

Our evaluation of this potential by numerically summing up the series is plotted in Fig. ^ in units 
of A/(2iV^). Just like the dilaton and contributions, the 2-form potential in Fig. || appears to be 
attractive. 

3.3.3 Asymptotics of the 2-form contribution 

As before the 2-form contribution V^^^ifl) only has poles along the imaginary-g axis. The plot of 
Vadj^'i) ^-loiig the positive Im g-axis is shown in Fig. 

The leading pole oiV^^{q) is found to be at mf = 5.1086±0.0001 (for A = 8x 10^ iterations), in 
agreement with [15]. The residue of this pole is i (0.051±0.003), such that the large-r T asymptotics 
of the 2-form contribution to the coordinate-space potential is 

2 A 0.26 ±0.02 

-T»i ~M 27rr 



-5.1086 TrrT 



(3.71) 
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Figure 9: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the 2-form field between the string world-sheets, plotted along the real-g axis (Im g = 0) in units of 
X/{2N^). The plot is for = 2 x 10^ iterations. 
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Figure 10: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the 2-form field between the string world-sheets, plotted along the imaginary-(7 axis (Re q = 0) in units 
of X/{2N^). The plot is for N = 10"^ iterations. 



Again this is an attractive potential. 



At small-r T/large-g the asymptotic solution of Eq. ( p.65| ) is 



which, when used in Eqs. ( |3.66D and ( p.2UD yields 



2^ ^ In 4^, (3.73) 



rT«i ^c'4vrr rT 



adj 

which is also attractive. 
3.4 The Graviton 

3.4.1 The potential and Einstein equations 

We want to find quadratic fiuctuations in the metric, qmn + ^mtv, induced in the presence of a 
straight string stretching from the UV boundary all the way down to the horizon. The most general 
form of the static metric with the 5*0 (3) rotational symmetry takes the form 

ds' = ^ [1 + h{r, z)] dr' + ^^^^ [dr' + r'dQj] + ^ jj^^ ^z' ' (3-74) 

It can be shown that there are consistent (non-singular) diffeomorphism transformations that can 
lead to the above metric. 

Since the coupling of the string world-sheet to graviton is 



Ana' 



^ j d'a ^j^'kMNdaX'^dbX'' (3.75) 

— * 

it follows that the 5-dimensional energy-momentum (EM) tensor of the string at X2 = in the bulk 



IS 



jMN ^J_Vl^ab ^(3) Q^x'I d,X^ (3.76) 

27ra' 

(see [41-43] for similar calculations). The non- vanishing components of the EM tensor are 

J'o = J: = ^ S^'\x) . (3.77) 

Therefore, in the presence of the string source, the Einstein equations can be written as 

= + 6^^ , (3.78) 

where Qff is the Einstein tensor, 6 is coming from the cosmological constant term, and = Att'^/N^. 

As usual we will work in units where Zh = 1 and -RAds = 1- Linearized Einstein equations ( p.78|) 
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for the metric ( 3.74|) are presented and simplified in Appendix ^ using momentum-space metric 
components 

h{r,z) = J ^e^'^-"f>(g,^) (3.79) 



for the rescaled fields 



A{r,z) = -^A, B{f,z) = -^B, h{r,z) = h. 3.80 

W W 

Kg 

As argued in Appendix ^ solving Einstein equations with additional gauge freedom left in the 
metric ( p.74|) allows one to set 

A + B + h = 0. (3.81) 

After lengthy but straightforward calculations presented in Appendix the equation for A (at the 
leading order in perturbations) can be brought to the following form: 

z{l- z^) (3 - 3z^ - z^) A. - [3 (1 - z^f - (3 + /)] A, 

-q^z [l-q^z' + 3z''] A = \q^z^ . (3.82) 

Again the boundary conditions are A{q, 2; = 0) = and A{q, z = 1) should be finite. 

Similar to the previous cases, the full graviton action evaluated on the classical solution of 
Eq. (^ is 

1 

Sk = ^J rfV V^r'hMNdaX^^d.X'^ = I ^ Air, z) , (3.83) 

(1) 

such that the graviton contribution to the adjoint heavy quark potential is 

1 

K5W = -^ / ^Air,z) (3.84) 
in coordinate space, translating into 



= / (3.85) 



in momentum space. 
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3.4.2 Solution of Einstein equations 

The series solution of Eq. (|3.82|) is 



A 



oo 
n=0 



,2 n+2 



(3.86) 



n=0 



with the recursion relations 

2 (2 + 2 n 



Pn 



3(2n + 3) (2n + 5) 

2(2n2-2n-l) 
~ 3(2n + 3) (2n + 5) 



q Pn-i - 



6-24n2 + g4 

3(2n + 3) (2n + 5) 
(2n-3) (2n-5) 
(2n + 3) (2n + 5) 



Pn-2 



Pn-i, 



135 

4n^-8n + l ^ 
^ 12n(n+l) ^ + 
13 - 16n + 4n2 
12n(ra + 1) 
So = 1, s-i = s_2 



P-3 

48- 



= 0, 

72n + 24n2-g^ 



(3.87a) 



■ q Sn-3 - ■ 
= S-3 = 0. 



12 n(n + 1) 
{n - 3) (n - 4) 



'?i-2 



n(n + 1] 



Neumann boundary conditions imposed on partial sums of the series ( |3.86|) yield 

J2n=0 Pn (2 n + 5) 



C^{q\N) 



En=o {2 n + 2) 



(3.87b) 



(3.88) 



Finally, the series representation for the potential is obtained by using Eq. ( p.86| ) in Eq. (|3.85|) , 
which results in 



V^Aq) 



X 



2 m 



E 

.n=0 



Pn 



2 (n + 2 



+ C^{q')Y. 



n=0 



2n+l 



(3.89) 



The numerical evaluation of the series is plotted in Fig. |TT|. The graviton naturally gives an attractive 
potential. 

The g-range of the plot in Fig. |ll] is limited by the convergence of the series in Eq. ( p.89|) 
both from above and from below, unlike the previous supergravity fields for which the convergence 
problems only appeared at large-g. At small (but non-zero) real q the numerical evaluation of the 
partial sums of the series in Eq. ( |3.89D appears to diverge, though the singularity moves towards 
g = as one increases the number of terms in the sums. One may then suspect a non-analyticity at 
g = in V^j{q). On the other hand one has to remember that in arriving at Eq. ( |3.82|) we have made 
a substitution (|A6|) , which is potentially singular at g = 0: the peculiar behavior of the numerics 
may be attributed to this potentially dangerous operation. To test whether this convergence issue 
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Figure 11: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the graviton field between the string world-sheets, plotted along the real-q axis (Im g = 0) in units of 
X/(2N^). The plot is for N = 10^ iterations. 

is an artifact of the numerics or a genuine singularity of V^j{q) at g = one may search for the 
solution of Eq. ( |3.82|) as a power-series in with the coefficient being some functions of z. A 
straightforward calculation yields 



A{q,z) 



(3.90) 



where, in order to fix the coefficient of the leading term one has to require finiteness of A{q, z 
at the order-g^. Substituting Eq. ( p.90|) into Eq. ( p.85| ) gives 



1) 



0) 



(3.91) 



We see that the graviton contribution to the heavy quark potential is finite at g = 0. Moreover, 
the series in Eq. ( p.90| ) can be easily continued, demonstrating that Vadji^l] is analytic at g = 0. 
Therefore the divergence described above is a numerical artifact and not a real physical divergence. 

3.4.3 Asymptotics of the graviton contribution 

The graviton contribution V^*^Jj(g) only has poles along the imaginary-g axis. The plot of V^*^|(g) 
along this axis is shown in Fig. |l^. 

The leading pole of Vadjifl) is at mf = 2.33 ± 0.01 (for N = 3000 iterations), also in agreement 
with [14,15], though with significantly lower precision, due to the complicated structure of the 
recurrence relations ( p.87|) for the series coefficients. The residue of the leading pole is i (0.37±0.01), 
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Figure 12: The contribution to the adjoint finite-T QQ potential in momentum space due to the exchange 
of the graviton field between the string world-sheets, plotted along the imaginary-q axis (Re g = 0) in units 
of X/{2N^). The plot is for N = 10^ iterations. 

giving the following large-rT asymptotics of the graviton contribution to the coordinate-space 
potential 



A 0.86 ±0.03 



-2.33 7rrr 



rT>l 



2 m 2 



TT r 



(3.92) 



This is also an attractive-potential contribution. 

At small-r T/large-g the asymptotic solution of Eq. (^.82|) is 



A\ ~ 

•^\gz»l 3^2' 



(3.93) 



which, when used in Eq. ( ^.85| ) yields 



Va%{r) 



■T<1 



2 A 1 1 



(3.94) 



which is, yet again, an attractive potential. 



3.5 The Net Result 

The net contribution of all supergravity fields to the QQ potential at the order- A/A^^ is given by 
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adding together Eqs. ( F^ , (|05D , j^M) , and ([ISSD , which yields 



2Ar2 



dz 



(^(g, 2;) + - t°(g, z)+Az^ Bo,{q, z) - ^(g, z) 



(3.95) 



The corresponding coordinate-space potential V^^ir^ 



(3.96) 



can be obtained from Eq. (|3.95|) using Eq. ( p.20| ). The total potential Vadjiq) is plotted in Fig. [13 
(lower line) in momentum space (along with the quark-quark potential to be discussed later in 
Sec. ^). Since all the four contributions to this potential were attractive, the resulting potential 
is indeed attractive. This is in contrast to the adjoint potential in perturbation theory, which 
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Figure 13: The net adjoint finite-T QQ potential (lower curve) and the net QQ potential (upper curve) 
in momentum space due to the exchanges of all supergravity fields between the string world-sheets, plotted 
along the real-g axis (Im g = 0) in units of A/ (2 N"^). 



was found to be repulsive (see e.g. [18]). As observed before, the adjoint potential (or the adjoint 
contribution to the unified total potential, as discussed in Sec. is of order- A/A*"^ both at weak 
and at strong couplings. We therefore conjecture that the difference between the strong- and weak- 
coupling regimes is in the sign of the interaction, making the repulsive weak-coupling potential 
attractive at strong coupling. 
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The large-rT asymptotics of the net potential V^^{r) is given by the mode with the hghtest 
screening mass, which, as follows from the above calculations (and from earlier works [14,15]) is the 
graviton, such that 

A 0.86 ±0.03 



adj 



rT>l 



-2.33 TrrT 



rT>l 



2 AT 2 



271 r 



(3.97) 



with the exponent in agreement with [11]. As we discussed above, if higher KK modes were allowed, 
than the large-rT scaling would be dominated by the field, and would be given by Eq. (|3.51|) . 
At small-rT, adding the contributions from Eqs. ( p.36|) , ( p.53|) , ( |3.73|) , and ( p. 941) yields 



v: 



adj 



rT<l 



-4 — In . 

Anr rT 



(3.98) 



Let us now return to the question of what happens to the potential V^^{r) as rT —t- 0. As 



adj 



follows from, say, Eq. (|3.95|) , the potential is obtained by integrating the contribution of various 
supergravity fields over the string world-sheet [dz/ z^). Our setup of static straight stings exchanging 
supergravity fields is valid only if the fields are weak. In coordinate space in the small-r T limit the 
fields are parametrically of the order 



A z 



(3.99) 



as follows from Fourier-transforming e.g. Eq. (|3.34|) (we include the coupling to the other string as 
well). The applicability region of the weak- field approximation is then defined by 



— - 

m r ^ ' 



which means that 



z <^ r ■ 
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X 



(3.100) 



(3.101) 



Therefore, assuming that corrections to our setup, coming presumably from string fluctuations, 
would regularize the divergence in the z integral at large- 2; (see e.g. [44] for similar phenomena), 
and putting z^ = I/^ttT) back into the expression explicitly, we see that the ^-integral in Eq. ( p.95|) 
is cut off by min {zh,r N^/ \} in the IR, such that 



V 



adj 
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m 47rr 
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again with the logarithmic accuracy. Therefore we expect 



Vafir) 
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(3.102) 



(3.103) 



which is finite in the rT — )■ limit. 
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4. The QQ potential 



We can use the results of the above calculation to find the potential between two heavy quarks 
immersed in the finite-T A/" = 4 SYM medium. While the decomposition into the two potentials 
corresponding to color group representations Nc and N"^ — Nc is possible along the lines of Eq. (2.6), 
as we will see momentarily both potentials obtained this way would correspond to the same string 
configuration in AdS space, and, therefore, would be equal to each other. Hence we define the QQ 
potential simply as 



— (TrL(O) TrL(f)), 



(4.1) 



in analogy to Eq. (|2.13| ). 





Figure 14: The quark-quark potential at finite-T at the oidei-X/N'^. 



To calculate V^^{r) using AdS/CFT we notice that now both strings attached to the quarks 
have to have the same orientation. Therefore the hanging string configuration from the left panel 
of Fig. H is now impossible, and the potential is given by the straight-strings configuration shown 
in Fig. |14|. Just as for the QQ system, the left panel in Fig. |1^ gives zero contribution to the 
potential after renormalization and subtraction of self-interactions. We are left with the right panel 
in Fig. Note that now both strings are oriented in the same direction, and can be parametrized 
as (cf. Eq. ( |3TTD ) 



Xi 



M 



(J 



X. 



M 



(r,0,2; 



a). 



(4.2) 



The effect of this difference is only felt by the 2- form field coupling to the string X2, which changes 
its sign, as follows from Eq. ( |3.54| ). Hence the action of the configuration in the right panel of 
Fig. 14 differs from that in Fig. 3 by the sign of the 2-form contribution. We therefore write for the 
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momentum-space QQ potential 



1 



A 



2 AT 2 



dz 



i^{q,z) + -i\q,z)-Az^Bo,{q,z)-A{q,z) 



(4.3) 



and for the coordinate-space potential 



VL{r) + V, 



The right-hand-sides of Eqs. ( [4 .31 ) and ( |4.4| ) contain terms given by the results of the above calcula- 
tions of the contributions to the QQ potential. The QQ potential from Eq. ( [4. 3D can be evaluated 
numerically using the above techniques and is plotted in Fig. ^ (upper line). It appears to be 
attractive, though somewhat weaker than the QQ potential. 

The large-rT behavior of V'^'^{r) is dominated by the graviton exchange, just like in the QQ 
case: 

A 0.86 ±0.03 



G 



(4.4) 
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(4.5) 



The small-rT limit is less straightforward to find. First of all, using Eqs. ( p. 361 ), ( |3.53| ), ( |3.73| ) 
and in Eq. yields 



(in the leading-logarithmic approximation) 



(4.6) 



-T<1 



This result can be interpreted as follows. The large-g asymptotics of the supergravity fields given 
by Eqs. (|3.34|) , (|3.52|) , (|3.72| ), and Eq. ( p.93| ), that we used in deriving the potential asymptotics in 
Eqs. (|33^) , (|3|), ( pT3D , and I^M) in Eq. (gj) comes from the part of the EOM for those fields 
that is independent of the curvature of AdSs. Hence this part of the string-string interactions is 
the same as in the flat space. In the flat space, however, we know that the interaction between the 
parallel strings oriented in the same direction should vanish since this would be a BPS system. In a 
way the zero of Eq. (|4.6|) provides a cross-check that we have properly accounted for all the relevant 
supergravity modes. 

The zero on Eq. ([4.6|) only implies that the leading-logarithm terms containing (1/r) ln(l/rT), 
which gave rise to the small-r T potential ( p.98|) in the QQ case, cancel among each other in the QQ 
case. This does not imply that the QQ potential actually goes to zero as r T — )• 0. To determine the 
correct rT — )■ asymptotics of the QQ potential we have to solve the EOM for the four relevant 
supergravity fields (0, Bmn, hjifN) at T = and integrate 
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2A^2 



dz 
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Putting Zh-i- ooin the dilaton EOM (|3T7D yields 



3 

- fz 
Z 



q (fi = z, 



with the solution 



(PT=o{q, z) 



TT 



2qz 



[L,{qz)-h{qz)]\. 



(4.8) 



(4.9) 



In arriving at Eq. ( [4.9|) we have required that ip{q^ 2; = 0) = and that if does not diverge more 
than linearly in 2; as z — )■ 00 (see Eq. ( |3.34| )). In Eq. ([4.9|), L2 is the modified Struve function and I2 
is the modified Bessel function. Due to the linear divergence at large-z, the contribution of (p from 
Eq. ( |4.9| ) to the QQ potential ( [4.7p is infinite: since we know that this infinity would cancel with 
contributions from other fields, we have to add those other contributions first before integrating the 
result over z. 

For the massive scalar t° the T = EOM in momentum-space is 



z. 



(4.10) 



Imposing the boundary condition t°(g,z = 0) = and requiring that t° is no more than linearly 
divergent in 2; as z —t- 00 (see Eq. ( p.52|) ) we obtain the solution of Eq. ( [4.10| ) 



35 
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TT 



+ —z^U{qz). 
z q 



(4.11) 



Here F is the generalized hypergeometric function. 

The zero-temperature EOM of motion for the 2-form field is 

z" El, + ^zBl,- (15 + q^ z^) Boz = z . 



(4.12) 



The series solution of this equation can be straightforwardly constructed (again imposing the same 
boundary conditions as for the scalars, with the exception that the 2-form field actually goes to 
zero as z ^ 00, as follows from Eq. (|3.72|) ). The summation of the series is more involved, but can 
be accomplished yielding 



qBl=''{q,z) = I,(q z) + -^ {271 [{8 + q'z')h{qz)-Aqz (l + 2 hiq z))] Y,{~tq. 



iqz 1 
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(4.13) 



where G is the generalized Meijer G-function. Despite several terms containing an i on the right- 
hand-side, Eq. ( [4.13|) gives a real-valued 2-form field. 
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Finally, the Einstein equations at T = 0, reduced to a single equation for A, read 



z{3-q'z')A..-3 (1 -q'z') A. - z [l - q' z'] A = ^q'z' 



(4.14) 



This equation can be solved by series expansion similar to how it was done in the non-zero tempera- 
ture case. However, the series can only be summed numerically. Instead we simply solve Eq. ( [4. 141 ) 
numerically, requiring that A{q, z = 0) = and that A maps onto the asymptotics of Eq. ( 3.93| ) at 
large-g z. 

Numerical evaluation of the resulting QQ potential at T = by the integration over z of the 
sum of all four contributions following Eq. (^4.7| ) yields 



Vt'S(9) 



-0T±0.1)^1. 



The corresponding coordinate-space potential is 
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(4.15) 



(4.16) 



We see that the net QQ potential at small-rT is attractive. We also observe that it is finite in 
the T — 7- limit, unlike the QQ potential of Eq. ( p.98| ), and does not require string-fluctuation 
corrections to remain finite in this limit. 



5. Summary and outlook 

We computed the adjoint quark-anti-quark potential, as well as the quark-quark potential in a 
SYM plasma. Both potentials are of order-A/A^^, and both are attractive. For comparison, notice 
that at weak coupling both potentials are also of oidei-X/N^. Moreover, the adjoint QQ potential 
is repulsive at weak coupling, and so is the the QQ potential with the quarks in the A*"^ — Nc 
representation. On the other hand, the QQ potential in the Nc representation is attractive at weak 
coupling. It appears that all the potentials, regardless of whether they are repulsive or attractive, 
become attractive at large 't Hooft coupling. 

To interpret this result we use the suggestion of [13] that the electric modes can be singled out 
non-perturbatively by identifying the CT-odd interaction channel. The quantum numbers of the 
supergravity fields contributing to the potentials calculated above can be found in [15] (see also [11] 
for a table of the quantum numbers of the relevant fields excluding t"). We see that the dilaton, 
graviton, and t° are CT-even, while Bq^ is CT-odd. We can thus identify the 2-form contribution as 
corresponding to chromo-electric modes in the gauge theory. The remaining dilaton, massive scalar 
and the graviton may correspond to either magnetic modes, or to even numbers of electric mode 
exchanges. The latter are probably not very important in the mix, since the chromo-electric (2- 
form) contribution by itself is smaller than the sum of all other contributions: hence it is likely that 
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the scalars and the dilaton correspond mainly to the chromo-magnetic modes on the boundary. The 
chromo-electric (2-form) contribution changes sign in going from QQ to QQ, which seems natural 
for the chromo-electric modes. In the perturbative weak-coupling limit, when the potentials are 
calculated at lowest order by two-gluon exchanges between Polyakov loops [18], the interaction is 
entirely due to chromo-electric modes. However, even at weak coupling it appears that the chromo- 
magnetic modes are less screened than the chromo-electric ones, though they do not couple directly 
to the heavy quark and anti-quark, and their contribution is suppressed by extra powers of the 
coupling coming from the loop diagrams needed to couple them to the quarks [13,45]. It seems 
that our above results suggest that at strong coupling the magnetic modes dominate over electric 
modes in the QQ and QQ potentials, making both of them attractive. The fact the the leading 
pole for the graviton is smaller than the leading pole for the 2-form field can be interpreted as 
magnetic modes being less screened even in the strong coupling limit. Small-coupling suppression 
of the magnetic contribution's coupling to the static heavy quarks by higher orders of the coupling 
is no longer an issue at strong coupling. Thus it seems natural that magnetic modes dominate at 
strong coupling (since they are screened less than electric modes and couple to quarks similarly) 
making both potentials attractive. 

For the adjoint QQ potential (and for the QQ potential in the A*"^ — Nc representation) we 
thus conjecture that in transition from weak to strong coupling the potential, while remaining of 
ordex-X/N'^, would change its sign. There should be no change of sign for the QQ potential in the 
Nc representation. 

In [11] the Debye screening mass was identified as the leading pole in the CT-odd channel. Of 
all the CT-odd supergravity particles, the one with the lowest location of the pole on imaginary 
axis is the axion [15], the EOM and, therefore, poles for which are the same as for the dilaton 
considered above. The authors of [15] have therefore identified uid = 3.4041 ttT as the Debye mass 
of A/" = 4 SYM plasma. However, it may seem a little peculiar that the axion does not couple to 
the string world-sheet at the leading (tree-level) order: how can electric modes in gauge theory not 
couple to a static heavy quark? It is possible (though seems a little hard to determine precisely) 
that the single axion can still couple to the string through some higher-order diagrams, with the 
axion-string coupling suppressed by additional powers of l/N^. If that was the case, the CT-odd 
part of the QQ potential at large-rT would be proportional to 

^ -5.W85nrT , 1 ^ \ „-3.404l7rr-T /r-i\ 

Indeed at very large distances, parametrically defined by rT ^ InA^^, the axion would dominate, 
with its mass being the correct Debye mass. It appears that more work is needed in order to 
determine the degree of suppression of a single-axion coupling to the string world-sheet and to 
eliminate the possibility of this coupling being zero. If the coupling is in fact zero, the leading pole 
of the 2-form field at 5.1085 ttT would be the correct Debye mass. We leave this for future work. 
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Other possible future improvements of our result may include performing similar calculations 
of the QQ and QQ potential in more QCD-like geometries [46-57]. Indeed the positions of the 
poles and the corresponding residues would be modified by the new geometries, along with other 
calculational details. One may even include finite quark masses (which are infinite in our present 
calculation) by inserting a probe D7-brane wrapped over 5*^ along the lines of [58] and having the 
strings dual to the quarks end on this D7 brane. After performing a more QCD-like calculations, 
one may hope to be able to perform a meaningful comparison of the results to the lattice QCD 
data [26,27]. 
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A. Solution of Einstein equations with a straight string as a source 

Linearized Einstein equations ( |3.78| ) for the metric ( p.74| ) read 

itt) SB + z{l + z^){?>A,-B,) + 2z{2 + z^)%-z^{l-z^)%^ + q^z'^\i = -'^z^ (Ala) 

(r r) SB + z[2{?> - z^) A. - {I ~ z^) + {I - z%, - z {I - z^) A., + z A 

-zrdl{A + B + \^)\ = (Alb) 

{9 9) 8 S + ;z [2 (3 - z^) A - (1 - z^) B, + {1 - z^) - z {1 - z^) A,, + z A 

-zdr{A+B + t})] = -z^ (Ale) 

3 

(zz) 8B + z[3{l + z'^)A,-2{2-z^)B, + {l + 5 z^) if, - 3 z {1 - z^) A,^ - z {1 - z^) \, , 

+ zq^B] = ~z^ (Aid) 

where we have omitted the (p 'P and r z components, which follow from the equations given above. 
As usual Az = dzA and Azz = dlA. 

Note that it appears difficult to perform Fourier transform on the r r and 9 9 components, since, 
for our 9- and ^-independent metric in spherical coordinates not all r-derivatives in those equations 
lead to powers of q in an obvious way. That is why we have, perhaps a bit sloppily, left some powers 
of r and dr in those two equations: these object should be understood now as some operators in 
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g-space. Subtracting r r and 6 equations from each other, and transforming the result back into 
coordinate space, we get 

-rd^,{A + B + h) + dr{A + B + h) = (A2) 

which imphes that 

A + B + h = pi{z)r^ + p2{z) (A3) 

with pi{z) and P2{z) some arbitrary functions of z. On the physical grounds one may require that 
the metric should fall off with increasing r at large-r: this would put pi and p2 to zero. Alternatively 
one can show that, employing the residual diffeomorphism invariance left in the metric ansatz ( p.74|) 
one can eliminate pi{z) and p2{z), obtaining Eq. ( |3.81| ), which in momentum space leads to 

A + B + ^ = 0. (A4) 

Using Eq. ( [A^ ) to eliminate 1) from Eqs. (|A1|) , and subtracting zz equation from tt, we write 



{tt)-{zz) -3(1- z^)A, + 3z{l- z^) A,, ~ zq^A-2zq^B = (A5a) 

(rr) 8B + z[{5-z^)A,-2{l-z^)B,-z{l-z^)A,, + zq^A] = ^z^, (A5b) 



which is now written entirely in momentum space. Solving Eq. ( [A5a| ) for B yields 



B = --^ [3il-z^)iA,-zA,,) + zq^A]. (A6) 

Substituting Eq. ( |A6| ) into Eq. ( |A5b| ) yields a 3rd-order differential equation for A. Solving this 
equation to express ^222 in terms of lower-order derivatives of A, we can use the result, along with 
Eqs. ( |A4| ) and (|A6|) , in, say, Eq. ( [Ala| ) to obtain Eq. ( p. 821 ) in the text (after some considerable 
algebra) . 
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